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O ■ Abstract 

O ; 

' We address the question of radiation emission from a perfect mirror that starts from rest and follows the trajectory 

1 z = — In(cosht) till i — > oo. We show that a correct derivation of the black body spectrum via the calculation of the 

' Bogolubov amplitudes requires consideration of the whole trajectory and not just of its asymptotic part. 

■ I. INTRODUCTION 

m 

' In a companion paper (Calogeracos 2002) hereafter referred to as I we addressed the question of emission of radiation 
from accelerated mirrors following prescribed relativistic asymptotically inertial trajectories. The fact that the mirror 
: moved at uniform velocity enabled us to define in and out states and employ standard time dependent perturbation 
theory. In the present paper we consider a perfect mirror starting from rest and accelerating for an infinite time along 
the trajectory 

z = g{t) = ln(cosh nt) (1) 

O , ' K 

I . The above problem has been considered extensively in the past, starting with the classic papers by Fulling and Davies 
^jpj' (1976) and Davies and Fulling (1977) (DF in what follows). The problem is of interest because of its (alleged) connec- 
• • . tion to radiation emitted from a collapsing black hole (Hawking 1975), (DcWitt 1975) and to the attendant thermal 
spectrum. DF have calculated the renormalized matrix element of the r„„ component of the energy momentum 
tensor, the latter being defined as 



X 



\2 



C3 ; Tuu = {ducf) 

(see e.g. Birrell and Davies (1982), equation (4.16)). The result is that asymptotically for t — > oo 



k2 



(T™) -> — (2) 

To derive (^ DF have used the asymptotic expression for the trajectory 

g{t) ^~t- Ae-2«* + B (3) 

where A, B are constants that can be readily determined from (|l|) . Equation (|^) shows that there is a constant energy 
flux at late times, which is interpreted (p. 249 of DF) as being analogous to the thermal energy flux found by Hawking 
(1975) in the case of a black hole. DF also calculated the Bogolubov amplitude (3{Lo,ijj') (and thus a{u),u)') as well) 
and the spectrum 

n{Lo)^ dLo'\l3{uo,uo')\^ e-'<'^+'^') (4) 
Jo 

(a is a convergence factor). They did find that the spectrum coincides with the black body spectrum, namely that 
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In this note we wish to re-examine the method used by DF to obtain the expressions for the Bogolubov amphtudes. 
Our findings are summarized in the concluding section, and there is no reason to anticipate them here. The motivation 
for the investigation is provided by the foUowing observations, (a) The Bogolubov a{LU,uj') and (3{uj,uj') amplitudes 
are by definition time independent quantities. Hence one should be sceptical as to the validity of using the limiting 
form (|^) ab initio without justification, (b) A careful calculation reveals that a certain term is missing from the 
starting equation (2.10a) of DF. Quite apart from the trivial nature of the error, it is not clear that the omission does 
not undermine the validity of the final result (^). 

As far as the physics of the problem is concerned one may well counter that in the context of mirrors acceleration 
continuing for an infinite time implies mathematical singularities and also entails physical pathologies associated, for 
example, with the infinite energy that has to be imparted to the mirror. In that sense the realistic problem is examined 
in I. In the present note we ignore such questions, take the premises of the early papers on the subject for granted, 
and concentrate on the calculation of the Bogolubov amplitudes and of the spectrum of the emitted radiation. In 
what follows we will show that asymptotic behaviour of i3{uj,llj') for large us' is 

Pico, J) « (c.')"^ + O ([oj')-"") {N > 1) (6) 

The Lo' integration in (Q) is then logarithmically divergent. The divergence signifies according to DF (p. 250) the 
production of particles at a finite rate for an infinite time. A similar divergence has been noticed by Hawking (1975) 
(p. 211) in the context of black holes and the interpretation offered is the same. One often refers to the ultraviolet 
divergence by saying that large uj' frequencies dominate. In that same asymptotic limit n(Lu) takes the familiar form 
of the black body spectrum (times the logarithmically divergent factor). Our analysis emphasizes two points: (a) the 
thermal result depends crucially on the fact that asymptotically the trajectory tends to a null line (the line i; = In 2 
in the case of (|^)), (b) for a mirror starting from rest and forever accelerating one must consider the whole trajectory 
and noi just the asymptotic portion of it. The truth of statement (a) is usually taken as common knowledge. However 
statement (a) is also taken erroneously to imply the negative of (b). 

Note: Notation and conventions follow that of I (and are largely in accordance with the papers cited in the 
Introduction). In section 2 we briefly go through the standard notation so that the paper be self contained. 

II. THE FORMALISM FOR A PERFECT MIRROR ACCELERATING FOREVER 

We introduce the usual light cone coordinates 

u = t — z,v = t + z (7) 

The massless Klein-Gordon equation reads 

dudv ^ 

Hence any function that depends only on u or w (or the sum of two such functions) is a solution of (^). Let 

z = 5(t) = -ln(cosht) (9) 
be the mirror's trajectory. Notice that for large t the trajectory assumes the asymptotic form 

g{t) w -t- e"^* + ln2 

Details are given in Appendix A of I; for the moment observe that the trajectory is asymptotic to the line w = In 2 
(figure 7 of I). In terms of the u, v coordinates the equation (|^) of the trajectory is written in the form 

u^fiv) (10) 

or alternatively 

V = piu) (11) 
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where the function / is the inverse of p. The construction of the functions / and p given the equation of the trajectory 
z = g{t) is explained in section 2 of I. In the particular case of @ the function / is given by 

f{v) = 1,« < 

/(w) = -ln(2-e''),0< u < ln2 (12) 



For future use we quote 



/'H = -j-^ (13) 



We take everything to exist to the right of the mirror. One set of modes satisfying and the boundary condition 

cjy{t,g{t))^Q (14) 

is given by (see I, equation (4)) 

1 

ipui(u, v) = — {e-x^i-iojv) - exp {-iujp{u))) (15) 

Another set of modes satisfying the boundary condition is immediately obtained from ( |l5| ) 

1 

(fuj{u,v) ^ — (exp {-iujf{v)) - exp {-iuju)) (16) 

2\/tTU! 



The modes ipuiu, v) of ( [Tsl) describe waves incident from the right as it is clear from the sign of the exponential in 
the first term; the second term represents the reflected part which has a rather complicated behaviour depending on 
the motion of the mirror. These modes constitute the in space and should obviously be absent before acceleration 

starts, 

a{uj') \Oin) = (17) 

Similarly the modes (pi^{u,v) describe waves travelling to the right (emitted by the mirror) as can be seen from the 
exponential of the second term. Correspondingly the first term is complicated. Regarding the out modes ( p^ ) the 
following point, which will be useful in the evaluation of the matrix element P{lu,uj') (|21|), ought to be noted. Recall 
the picturesque form of such a mode (figure 6 of I), where for a perfect mirror Rj^ — 1, Tfj — 0. Clearly in the 
coordinate range u < 0,v < the mirror is stationary and thus the modes trivialize. During the subsequent motion 



of the mirror the variable v in the back-scattered wave in (16) is restricted between and In 2. These remarks will 



reflect on the range of integration in ( 23 ) 



The (puj{u, v) modes define the out space and 

a{uj) \Oout) = (18) 

The state \Oout) corresponds to the state where nothing is produced by the mirror. The two representations are 
connected by the Bogolubov transformation 

/•OC 

a{uj)= duj{a{Lu,Lu')a{uj')+(3*{uj,uj')a\uj')) (19) 
Jo 

Using ( p9| ) and its hermitean conjugate we may immediately verify that the expectation value of the number of 
excitations of the mode (w) in the \Oin) vacuum is given by 

/>oo 

{Oin\ N (uj) \Oin) = duj' \/3{uj,uj')\^ (20) 







The matrix element P{uj,uj') is given by (Birrell and Davies 1982), equations (2.9), (3.36) 

P{iu,Lu')=-i J dzip^,{z,0)-^^^^{z,0) + i J dz(^-^^^^,{z,0)^^^{z,0) (21) 
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The integration in ( pl| ) can be over any spacelike hypersurface. Since the mirror is at rest for t < the choice t = 
for the hypersurface is convenient. The in modes evaluated at i = are given by the simple expression 



1 

LPuj{u, v) — — r= (exp(— icju) — exp {—iuju)) (22) 
2\/t:uj 



I.e. expression (|l|) with p{ u) =: u (corresponding to zero velocity). The (p modes are given by ( |l6| ) with / depending 
on the trajectory. Relation (|2l|) is rewritten in the form (the endpoints of integration shall be stated presently) 



(3{cu,cu') = -i dz^^\e-"^'' ~e'^''}-^{f'e-'^^ -e"^'}+ (23) 



2v7rw 

dz^^= (e"'"'^ - e'"'4 je"*"^ - e'^^ 



Bearing in mind the comments preceding (|l^) regarding the range of integration, the above expression may be 
rearranged in the form 

1 /"'"^ r ~i (lli — u}'] r°° r 1 

(3(uj, uj') = ^= / dzle'"^' - e-"^ ^ I { w'e"^"/ - uofe'"^^ } + ^ / i e*" ^ - e"'"'^ i e'"^ (24) 

47rVLi^ JO ^ ^ 47rvaJa/ Jo ^ 

where the limits of integration are displayed. The first and second integral in the above relation will be denoted 
by (3i{uj,uj') and Piii{uj^lo') respectively (in accordance with the notation in I). Notice that the second integration 
is of kinematic origin and independent of the form of the trajectory. We also quote the expression for the a{uj,uj') 
amplitude 

a{uj,uj')^ij dzip^,{z,0)—(p*^{z,0)-ij dz \^—(p^- {z,0) J ipl{z,0) (25) 
Observe the unitarity condition (see I for details) 

duj (a (wi, w) a* {lj2, uj) — f3 (wi, w) f3* (w2, ui)) = S {oji — L02) (26) 



which is a consequence of the fact that the set of in states is complete. Notice that identity (69) of I is not valid 
because the set of out states is not complete (in the present case where the trajectory has a v asymptote). On the 
question of completeness see the remarks preceding equation (35) of I. 

Recall from I that we can introduce quantities A{lj,lj'), B{u!,u!') that are analytic functions of the frequencies via 

A(lj,lj') ^, B(ijj,uj') 



The quantity B{uj,uj') is read off ( |2^ ) (and A{uj,lo') from the corresponding expression for Q!(a;, cj')). From the 
definitions of the Bogolubov coefficients, the explicit form (^3|) of the overlap integral and expressions p5|) and ( [l^ ) 
for the field modes one can deduce that 

B*{uj,uj') = A{^iu,uj'),A*{uj,uj') = B{~u,u') (28) 

The above relations allow the calculation of a(Lu,Lu') once /3(w,w') is determined. 



III. COMPARISON BETWEEN REALISTIC TRAJECTORIES AND TRAJECTORIES ACCELERATING 

FOREVER 

A mirror trajectory involving acceleration forever cannot be considered as realistic in contrast to an asymptotically 
inertial trajectory. The former type of trajectory is of interest in that it can be taken as a simplified analog of black 
hole collapse. The two trajectories are radically different from a physical point of view and this is reflected in the form 
of their spectra. In I we considered a mirror starting from rest, accelerating along the trajectory (|^) till a spacetime 
point P (labelled by a coordinate w = r < ln2 in {u,v) coordinates), and then continuing at uniform velocity. We 
showed that in that case the Bogolubov amplitude squared |/3(a;, w')] behaves asymptotically for large uj' as {uj')~ ■ 
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In this work we set out to show that if acceleration is to continue forever then goes as 1/w' (see (|)). It 

is thus not the case that the second problem can be simply considered as the r ^ ln2 limit of the first. Although 
this may not come as a surprise from a physical point of view, it is certainly interesting to see where the dichotomy 
occurs mathematically. To this end we return to the problem examined in I where the u') amplitude is given by 
( p^ ) above with In 2 replaced by oo (cf also equation (34) of I), with the understanding that the function f{z) stands 
for facc{z) given by ( |l2[) above when z < r, and for /o(z) given by (13) of I when z > r (in the last equation B stands 
for the velocity (3p at point P where acceleration stops and the constant C is adjusted so that velocity be continuous 
at P; see equation (49) of I) . We introduce the Fourier integrals 

POO 

h= dze"^''-°''e-''^f^''> (29) 
Jo 



dze" 



(30) 



dze 



'iuj z — az —iLuf(z) 



dze-"^''-"'f'iz)e-''^f^'''> 

Jq 

where a is a small convergence factor taken to zero at the end of the calculations. Thus the /-dependent part of 
amplitude (35) of I is a linear combination of the above four integrals. The treatment that follows is equivalent to 
that of section 3. A of I, however it makes clearer the contrast between the two cases (realistic trajectory vs one that 
accelerates forever). 

To examine the asymptotic behaviour for large uj' of the above integrals we use simple integration by parts (see 
Bender and Orszag 1978), p. 278 according to the formula 

''If'' 

/ dzF'(z)e*"^ (31) 

a J a 

which works provided the quantities appearing in the right hand side are well defined. To treat Ii of ( p9| ) we split the 
integral to + and apply ( pl| ) twice. The endpoint contributions from infinity vanish due to the convergence 
factor. The endpoint contributions at z = r cancel out in pairs: fo{r) cancels with facc{r) (both accompanied by a 
factor 1/lu') and /o(r) cancels with facd''') (both accompanied by a factor {uj')~ ). The origin of these cancellations 
is a direct consequence of the nature of the trajectory considered. The first cancellation reflects the fact that the 
trajectory itself u = f{v) is continuous, and the second reflects the continuity of the velocity at P (see equation (15) 
of I for the general connection between f'{v) and velocity). Similarly endpoint contributions at z = cancel with 
corresponding terms originating from a large uj' expansion of the second (trajectory independent) term in ( p^ on the 
same continuity grounds as above. Thus integral Ii goes as {u') ^ and observing the prefactors in (|2j) we deduce 

that its contribution to f3{uj,uy) goes as (w')"^^^. 

The integral I2 given by (pOj ) can be handled in a similar way. However because of the presence of f'{z) we can 
integrate by parts only once if the endpoint contributions are to cancel out (recall that /^{r) — facd''') would require 
continuity of the acceleration). Thus I2 goes as (w')^ • Because of the corresponding prefactor in (p3) is one power 
of Lo' smaller, the final contribution to the amplitude is again of the order {uj') 

We return to ( p^ with the upper limit of the first integral now being equal to In 2. It is manifestly obvious that 
there is now no room for the subtle cancellations exhibited above. We can still try integration by parts to see whether 
we can possibly arrive at any conclusions. To this end consider (for example) I2 (with the upper limit now equal to 
In 2). Recall that according to ( p^ ) 

and according to ( p^ ) 

2 — e^ 



dzF{z)e" 



F{z) 



lUJ' 
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Then 



/•In 2 , z 

Jo (2-6-)^+^" 

If we now attempt to apply property (|3l| ) the upper endpoint contribution diverges. Hence we are unable to extract 
a further power of l/uj'. It is clear that to obtain an asymptotic estimate in the present problem we have to resort to 
other methods. 



IV. CALCULATION OF THE BOGOLUBOV AMPLITUDES 



The strategy we adopt in handling (24) is as follows. The first integral will be evaluated via an asymptotic expansion 
in negative powers of uj' , which will in fact show that the a;' integration in (^) is logarithmically divergent. The second 
integral (3iii{uj,uj') in (|4|) is readily evaluated (see also I) 

(iiii{LO,Lo')^ . V^ (^-^OC(^ + c^O (32) 

where the function C and its complex conjugate ^*are defined in Heitler (1954), pages 66-71: 

C(x) = -i I e'^^c/K = P- - ni5{x) (33) 



As explained in I it is only the first term in (|3^) that is operative as far as the calculation of the /3(aj,aj') goes. 
(The 8 proportional term is only relevant in the calculation of the a(a;,a;') amplitude via relations (|27|), (p^).) Thus 
asymptotically in the large limit 

/3„7(^,^')«— ^ (34) 

We turn to the first integral I3j{lo,uj') in (p^). Rather than dealing with four integrals we perform an integration 
by parts to get (E4) in the form 



/3/(w,tj') = \ — / dze-"^^'^^^-"^''' + ^sin(tjln2)e-^"^('"2) 

27r V UJ Jq 27r ^/Zj^ 



(35) 



This is the same integration by parts that was used to obtain (54) of I and also the one that is used in DF to go 
from their (2.10a) to (2.10b). The second term in ( ^5|) oscillates rapidly since the exponent tends to infinity. Thus 
the term tends distributionally to zero and it may be neglected as in DF (also it is one power of to' down compared 
to the first term). Of course this term was kept in (47) of I since it gives a finite contribution for r < ln2 (recall that 
&t V = r the trajectory considered in I reverts to uniform velocity). So asymptotically we are entitled to write 

1 I T /-In 2 1 I J p\n2 



P,[uj,uj')^-^J^ / d.e— /(^)— ^--^V- / dze— (2-e-y'" (36) 



2tt y UJ Jq 2tt y uj jq 

To bring the singularity to zero we make the change of variable 

z==ln2-p (37) 

and rewrite Pj(uj,uj') in the form 



We isolate the integral 



In 2 







' dpe'^'p {l-e-py" (39) 
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To obtain the asymptotic behaviour of ( |39D we adopt the standard technique of deforming the integration path to 
a contour in the complex plane (see (Bender and Orszag 1978), chapter 6); see also (Morse and Feshbach 1953), p. 
610 where a very similar contour is used in the study of the asymptotic expansion of the confluent hypergeometric. 
The deformed contour runs from up the imaginary axis till iT (we eventually take T — *■ cxj), then parallel to the real 
axis from iT to ir + ln2, and then down again parallel to the imaginary axis from iT + \n2 to In 2. The contribution 
of the segment parallel to the real axis vanishes exponentially in the limit T oo. We thus get 

1 = 1 dse-^'' (1 - e-'^y^ - i / dse'"'(i"2+,.) _ g-in2-zsy- ^ 
Jo Jo 



= ij dse-'^'' {l-e-'^Y"' -iT"^' J dse"'^'" (^1 - ^ j (40) 
In the limit of cu' large the main contribution to (E^ comes from the s « region. We expand e** in powers of s. We 



thus approximate one factor of the first integrand in (40) by 

1 - e~'' « is, {tsf^ = e-'^'^^h''^ 

where we set 

,took the branch cut of the function p"^ to run from zero along the negative x axis, wrote p"^ = exp {iuj (Inp + i2NTi)) 
and chose the branch iV = 0. For the second integrand we get 

1 is 



Thus 



/ = ie-""^'^ / dse-'^'^s"^ - i2*('^'-'^) / dse'^'' = (42) 







le 



V{l + iuj) i2*('^'~"') 



We substitute (H) in (||) and get 

/3z(^, J) « (c^')""" e'^-'^T (1 + ^c.) + (43) 

Observe the crucial cancellation of the second term in ( |4^ ) with (|3^), the end result being 

As already mentioned at the end of the previous section, having determined /3(a;, oj') allows one to determine ot(uj, oj'). 
This is the typical form of the /3(ci;,w') amplitude that leads to the thermal spectrum; see e.g. DF. Indeed by taking 



the modulus of (44), squaring and using the property 

|r (1 + iy)f' = T^yl sinh (Try) 

we get the black body spectrum (0). 
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V. CONCLUSION 



The objective of the paper was to prove that the Bogolubov amphtude /3(w,cj') has the asymptotic form and 
that the radiation emitted has the spectrum of a black body. Before enlarging on the conclusions we should elucidate 
one technical point in connection to the classic paper by Davies and Fulling (1977). It seems that the term (3jij{uj, u)') 
is unaccountably missing from (2.10a) of DF. Its existence is necessary if the unitarity conditions are to be satisfied. 
The need for the (3iii term may be seen in a trivial example based on the results of I. Consider the limiting case 
of a mirror perpetually at rest. Then Pi{u;,uj') ((52) of I) disappears (formally r = 0) and Pni is instrumental in 
cancelling (3ii of (53) evaluated at r = 0, /3p = so that the total emission amplitude /3 vanishes (as it has to). The 
omission has been pointed out by Walker (1985) who however did not pursue the matter any further. On the other 
hand the present derivation shows that its existence is crucial in cancelling the non-thermal part of the Pi{lo, lo') term. 

There are various arguments in the literature in favour of the black body spectrum in the case of trajectory (|^). 
In the present note we are concentrating on a proof based on the Bogolubov coefficients. These quantities are by 
definition time- independent, and in this context the question as to where and when the photons are produced simply 
does not arise. In the same vein it is totally arbitrary to assert from the start that one specific part of the trajectory 
(in the present case the one near the asymptote) is more important than other parts. It is certainly true that were 
it not for the singularity on the v asymptote the thermal spectrum would not arise. However the main conclusion of 
this note is that the correct derivation of the thermal result requires the consideration of the complete trajectory and 
not just of its asymptotic part. In technical terms the function f{z) in ( ^ ) cannot be approximated by its asymptotic 
expression. These remarks are strengthened by (i) the aforementioned role of the missing term (which appears as if 
it originates at t = 0), (ii) the failure of the short distance expansion from the v asymptote as demonstrated in the 
Appendix. All this is in accordance with one's quantum mechanical intuition. One's classical instincts might dictate 
that roughly speaking the small amount of time spent near the origin would have an insignificant effect compared to 
the infinitely long time spent near the asymptote. However such loose statements are misleading in connection to the 
quantum mechanical calculation of global [time independent) quantities. Similarly attempts to distinguish between 
" transient" and " steady state" radiation at the level of the a and /3 amplitudes are bound to fail; the emphasis in 
the literature on the importance of the asymptotic part of the trajectory has unfortunately led to such statements. 
Having said all that, one can certainly enquire about the matrix elements of local field quantities as it is indeed done 
in DF in a most illuminating way. This is however quite distinct from the calculation of Bogolubov amplitudes. 

Mention must be made to the work of Carlitz and Willey (1987) where a mirror accelerating from the infinite past 
to the infinite future is considered. The authors do get the black body spectrum, their amplitudes do satisfy the 
Bogolubov identities, and quite clearly the time i = plays no special role in their problem which is quite different 
from ours and our comments do not apply to their work. 
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Appendix: Short distance expansion from the horizon and large frequencies 

Let us consider the integral ( p9[ ) and expand in small p. To first order 

An 2 

= / dpp^^e"^'" (45) 
Jo 

The integral (^5|) can be performed exactly in terms of the confluent hypergeometric and the asymptotic limit of large 
uj' may be examined afterwards. Let us make the change of variable p = t In 2 in (Kq) and rewrite 



/(I) - (In 2)*"+^ f\ 
Jo 



dte"^''''H"^ = (46) 

= (In 2)*"+^ — - — M (1 + iuj,2 + iuj, iJ In 2) 
iuj + \ 

(where M is the confluent hypergeometric function). An alternative way of presenting the above result is in terms 
of the incomplete gamma function 7* by exploiting the connection of the latter with the confluent (see e.g. Tricomi 
(1954)) 
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We can now examine the asymptotic limit of (46) for large oj' . The asymptotic limit of the confluent M{a,b,i \z\) for 
large values of \z\ is given by item 13.5.1 of (Abramowitz and Stegun 1972) (z = iuj'\n2). In the case b = a+1 some 
simplifications occur and we get 

M {l + iu,2 + iuj,i\z\) ^ -{1 + tuj) e'l^l -^ + ir{2 + iu) (48) 

\z\ 

(other terms are down by higher powers of 1/ \ z\). The second term of the above relation combined with the prefactors 
of ( ^6| ) does feature the F (1 + iu) factor characteristic of the black body spectrum. However the presence of the 
first term spoils the thermal result (notice that both terms are of the same order in u;'). In other words the result ( |4^ ) 
is the correct answer to (^), which in turn is the wrong approximation of the original amplitude. The error is hardly 
surprising, since the term has been omitted and the term /?/ has been approximated in a totally unsystematic 
way. 

To demonstrate that expansion in powers of p is a non-starter we go one step further in the expansion of ( |39| ) 
retaining terms of order . Denoting this second order approximation to / by J*-^^ we get 

1^^) = -— r^pp'^+^e'^'p (49) 
2 Jo 

The calculation of ( ^ ) proceeds along exactly the same lines as that of ( ^6| ) and yields 

/(2) = _!^(ln2)*'"+'^— M(ia; + 2,ia; + 3,iw'ln2) (50) 
2 iLu + 2 



In the — > oo limit 



i\z\ 



M{ioj + 2, zc^ + 3, z \z\) « -z (2 + ZL^) ^ - r (3 + luj) (51) 

\z\ \z\ 



The second term in ( |5l| ) ought to be neglected compared to the second term in (|48|) since the former is down by one 
power of Lo' . The first term in ( pT| ) is of the same order in uj' as the two terms in ([48|) and thus ought to be retained. 
The pattern persists to all orders, and is due to the structure M{iuj + n,iuj + n+l, iuj' In 2) {n integer) of the confluent 
in the present problem (the second argument being equal to the first plus one). Thus an expansion in powers of p 
gives contributions of the same order of magnitude (i.e. O (1/w') and hence is of little use. 

In retrospect it may seem quite surprising how one obtains the correct result starting with the wrong expression 
(p5|). Let us return to that expression which (as detailed above) is the first order approximation in a short distance 
expansion from the asymptote. The integral is often handled as follows (see for example (Birrell and Davies 1982) p. 
108. One rescales the variable and writes the integral in the form 

/(I) = {cj'y'^ r dpp'^e'P (52) 



One now simply sets Ct''ln2 = oo, changes variable p = ia and rotates in the complex plane to get ( p2[ ) in the form 

poo 

/(I) ^ (cj')"'" / dfre-V'" (53) 



"'0 

Notice that setting ti;'ln2 = oo certainly does not amount to a systematic expansion in (ut') ^. The a integration 
yields r(l + iuj) and one thus obtains the form for the (3 amplitude leading to the black body spectrum. On the other 
hand the integral (^) can be performed exactly (see (^6[) above) and the asymptotic estimate for large lu' (recall that 
we are chasing the ultraviolet divergence) is given by (p:8[) . The reason for the discrepancy lies in the fact that one 
should first evaluate the integral in terms of the confluent and then take the lo' oo limit rather than take the limit 
first. This rotation in the complex plane stumbles upon the Stokes phenomenon for the confluent (different limits 
for 1 2; I ^ 00 depending on argz)). In other words it appears that Birrell and Davies have made two self-cancelling 
mistakes (wrong approximation to the amplitude and wrong evaluation of the integral). 
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